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WO 98/31142 PCT/US97/24136 

DEFECT CHANNEL NULLING 



RELATED APPLICATION 

5 This application claims the benefit of U.S. Provisional Application No. 

60/035,763, filed January 6, 1997. 

FIELD OF THE INVENTION 

This invention relates to electronic image enhancement and recovery, and more 
10 particularly to a method and apparatus which compensate for the effects of defects of the 
media on which the imagers stored. . 

BACKGROUND OF THE INVENTION 

Ever since the first image of an object was captured on film, a serious problem 
15 became apparent: imperfections in and on the recording medium itself which distort and 
obscure the original image that was sought to be captured. This non-image imperfection 
problem continues to plague the field of image capture and reproduction to the present day. 
These imperfections occur in many forms including dust, scratches, fingerprints, smudges, 
and the like. 

20 The film industry has been concerned that the problem caused by these non-image 

imperfections may jeopardize the long term future of analog images. Notwithstanding the 
significant efforts that have been made to solve the problem, it nevertheless persists. This 
is particularly true with respect to enlargements and high resolution scans. Thus, the 
problem is becoming even more acute and recognized as resolutions increase. 
25 Furthermore, multimedia applications have brought widespread attention to the problem 
with the increase in film scanning for computer applications. 

The non-image imperfection problems occur more frequently with negatives than 
with transparencies because (1) viewing the negative requires a light to dark ratio gain (also 
known as a "gamma") of .about two to one; whereas viewing a transparency is less 
30 demanding and (2) filmstrips are subject to more contact than are mounted transparencies. 
Some imperfections may be present in fresh photographic film, for example, media surface 
waves or microbubbles in the emulsion. Improper processing may also introduce 
non-image imperfection problems, for example mild reticulation or residual unbleached 
silver. Some non-image imperfection problems are introduced through improper handling 
35 of the image media which may leave oily fingerprints and other surface contaminants or 
cause physical damage such as scratches which distort the view of the image. Improper 
storage may also leave the film further exposed to defect inducing environmental influences 
such as dust and chemicals. Archival polypropylene sleeves employed to protect negatives 
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contribute to the problem by leaving hairline surface scratches as the negatives are pulled 
out and replaced into the sleeves. As illustrated by the above examples, non-image 
imperfections may arise out of one or more of the following: film emulsion defects; 
improper processing; and improper handling or storage which may cause extraneous matter 

5 accumulation or scratches on the image media. 

Numerous measures have been developed in the art in an attempt to address these 
problems, particularly with respect to captured high resolution images. Generally 
speaking, these measures are in the nature of prevention and labor intensive minimization of 
non-image defects. One approach has been that of preventing introduction of non-image 

10 defects in the development process itself. Expensive anti-static equipment, including spray 
machines to neutralize dust attracting charges, are employed by some photo finishers. 
Photo finishers also attempted to solve some of the problems by employing diffuse light 
source enlargers that help reduce the effects of refraction of the light due to imperfections in 
the photo finishing processes. 

15 Another approach used in the an involves minimizing the effects of these 

imperfections once they are present by various correction techniques, most of which are 
manual and thus highly labor-intensive and expensive. For example, during the photo 
finishing process, a highly trained individual might spend a great deal of time using various 
spotting dyes and an extremely small spotting brush in an effort to essentially paint out the 

20 imperfections in a reproduced or printed image. Another technique is to wipe on or 
immerse the negatives in a light oil in an attempt to optically fill scratches. 

As the art has developed, various attempts have been made to automate the 
correction process, particularly with respect to digital image systems. In such systems, 
once an imperfection has been detected, various "fill" algorithms are used to correct the 

25 image at the imperfection site. Nevertheless, heuristics or human intervention has been 
required to detect the imperfections with a subjective threshold. Typically the area 
identified to be corrected in this manner is much larger than necessary, in part due to 
application of these subjective criteria for detecting defective areas. 

Automated methods have also been developed for detecting imperfect areas in 

30 recording media, as described in German patent 2821868.0 published Nov. 22, 1979, and 
entided "Method and Device for Detecting Recording and Counting of Mechanical Damage 
to Moving Bands, for Example Films." The approach discussed focuses on determining 
the quantity of defects and shutting down processing if that measured quantity exceeds 
some predetermined maximum level of defects. In this system a source of infrared energy 

35 impinges upon the film medium. A scanned infrared image of the film in question is then 
taken by sensors detecting reflection of the infrared energy from the film surface. 
However, several limitations are present in the system disclosed by the patent 
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First, its purpose is not to correct the effects of such detected film defects present on 
the film image. Instead, the system is implemented simply to monitor the prevalence of 
these defects in an automated photographic development process whereby the process can 
be automatically shut down if the defect rate exceeds a prescribed level. The optical 
5 infrared path is a reflective one from the infrared source to the infrared sensor which is 
typically different from the other sensors utilized for image processing. The infrared image 
is not recorded in registry with any other scanned images from the visual portion of the 
electromagnetic spectrum. Registry refers to the precise alignment of two related images, 
plates, or impressions so that they are held in position relative to each other. The fact that 
10 this prior art system does not record the infrared image in registry with any other images 
taken of the film is a disadvantage which in and of itself renders it extremely difficult to 
subtract out the effect of such imperfections noted in the infrared image from similar defects 
present in the visual record of the image. 

In another prior system disclosed by the present inventor, a method which 
15 compensates for the effects of storage media defects on image data is disclosed by deriving 
from the medium separate images in the red, green, blue and infrared portions of the 
electromagnetic spectrum, each corresponding to the image stored therein. As disclosed in 
U.S. Patent No. 5,266,805 issued to the present inventor on Nov. 30, 1993, red, green, 
blue and infrared light is sequentially directed at one side of the film by means of a light 
20 source and color filter wheel Corresponding red, green, blue, and infrared images formed 
by that portion of the light being transmitted through the film are digitally captured from the 
opposite side of the film. The images are preferably captured in registry to facilitate 
subtracting out the effect of imperfections at locations in the infrared record from 
corresponding locations in the red, green, and blue images. The imperfections may either 
25 substantially reduce or totally occlude the infrared light. However, remaining portions of 
the medium having the desired image without such imperfections are essentially uniformly 
transmissive to infrared light. These imperfection-free portions have variable transmissivity 
in the visual spectrum which is determined by the image developed on the film. 
Accordingly, the infrared image may serve as an indicator or map of the spatial position of 
30 these non-image imperfections on and in the medium, thereby allowing determination of the 
location and removal of the defects so that the underlying desired image may be recovered. 

In order to remove film defects from a scanned film image, an infrared image record 
is used as a norming channel by dividing each pixel in a visible channel (red, blue, or 
green) by the corresponding pixel in the associated infrared channel. Although this 
35 improves the image, variations in scanner contrast across the density range of the film, 
variations in image focus, variations in illumination aperture and variations in the way 
defects respond to infrared light cause imbalances between the visible and infrared records 
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leaving defects in the resulting image after the division. This process typically decreases 
the level of defects caused by storage of the image on a medium to within approximately 
10% of their complete elimination. Although greatly improved relative to the uncorrected 
original image, the nulling of defects achieved via division with the infrared image record is 
5 not sufficient for some purposes. It is, therefore, desirable to have a method insuring that 
defect nulling approaches 100% elimination for these purposes. 

When a transparent film is viewed, several factors attenuate light. The actual image 
attenuates light by dye absorption. In color film, each dye is specific to one color (red, 
green, or blue), and none of the three color dyes attenuate infrared light. Defects and 

10 imperfections attenuate incident light by bending the light outside the angle seen by the lens. 
For example, scratches, fingerprints, small flecks of dust, air bubbles, and other 
imperfections refract, or bend, light away from the lens, and thus appear as darker areas in 
the stored image. This refraction by defects is nearly the same in infrared light as in visible 
light; therefore, an infrared image shows the location and magnitude of defects without the 

15 visible image. For any given pixel, if the visible measurement is divided by the infrared 
measurement, any attenuation caused by defects will be divided out, in theory leaving only 
the attenuation from the visible dye image. The theory of infrared division is surprisingly 
powerful. 

In practice, however, it was found that defects did not attenuate infrared light 
20 exactly the same as visible light. The discrepancies between theory and practice are 
attributable to several factdrs. In particular, the index of refraction is slightly different 
between infrared and visible light. For example, (1) some dust flecks, even though very 
thin, may absorb some light; (2) the infrared and visible light sources in the scanner may 
not be precisely aligned causing rays bent by a defect to be "seen" somewhat differently 
25 under infrared light than by visible light; and (3) the imaging system in the scanner typically 
has less sharpness when using infrared as compared to visible light. Because of these 
discrepancies, if the visible record is simply divided by the infrared record as in the method 
disclosed by above-referenced U.S. Patent 5,266,805, defects are greatly attenuated but not 
entirely eliminated. That method leaves a ghost, or residue, of the defect that still requires 
30 some type of manual intervention to fix. It is an object of the present invention to eliminate 
manual intervention in correction of defects. 

Accordingly, it is an advantage of the present invention to enhance image recovery 
and eliminate effects of media surface defects on images stored in that media. 

To achieve this and other objects which will become readily apparent upon reading 
35 the attached disclosure and appended claims, an improved method of and apparatus for 
defect channel nulling which significantly reduces the effect of defects on the image 
recording medium is provided. Additional objects, advantages, and novel features of the 
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invention will be set forth in part in the description which follows, and in part will become 
apparent to those skilled in the art upon examination of the following, or may be learned by 
practice of the invention. The objects and advantages of the invention may be realized and 
attained by means of the instrumentalities and combinations particularly pointed out in the 
5 appended claims. 

SUMMARY OF THE INVENTION 

A method and apparatus for nulling the effect of media defects introduced by storing 
an image on a medium. Red, blue, and green image signals are collected by passing red, 

10 blue, and green light through a medium containing an image and having some surface 
defects. The collected red, blue and green image signals are affected by the presence of 
image medium defects. A defect signal indicating the defects of the image storage medium 
is received from the same defective image storage medium. The defect signal is in register 
with the image signal. Gain is applied to the defect signal to generate a normalized defect 

15 signal. The normalized defect signal is then subtracted from the image signal to generate a 
recovered image signal. The defect signal and recovered image signal are compared, and 
from this comparison a second gain is derived which is proportional to the degree that the 
deterioration indicated by the defect signal is included in the image signal. This second gain 
is then applied to the defect signal and this" now normalized defect signal is then subtracted 

20 from the image signal to generate a recovered image signal to minimize the effects of defects 
due to storage, handling, processing or which are inherent to the medium prior to recording 
of any image. 

BRIEF DESCRIPTION OF THE DRAWINGS 

25 Fig. 1 is a schematic representation of a representative system for implementing the 

invention. 

Fig. 2 is a flow diagram indicating the overall processing sequence for defect 
channel nulling. 

Fig. 3 is a flow diagram indicating the sequence of variable initialization. 
30 Fig. 4 is a flow diagram of the pyramid building process. 

Figs. 5A, 5B, 5C and 5D are diagrams which illustrate the process of pyramid 
downsizing. 

Fig. 6 is a diagram which illustrates the process of pyramid upsizing. 
Figs. 7A, 7B, 7C, and 7D are diagrams which illustrate the process of pyramid 
35 upsizing. 

Fig. 8 is a diagram depicting the process of differencing image-related signals. 
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Figs. 9A and 9B are flow diagrams indicating the steps in the process of correlating 
the image and defect pyramids. 

Figs. 10A and 10B are flow diagrams indicating the steps in the process of 
correlation histogram generation. 
5 Fig. 1 1 is a graph of the correlation histogram. 

Figs. 12A and 12B are flow diagrams indicating the steps in the process of 
calculating the correlation cutoffs. 

Fig. 13 is a flow diagram detailing the steps in the process of calculating the cross 
correlation and autocorrelation. 
10 Fig. 14 is a flow diagram indicating the steps in the process of calculating the 

clamped element. 

Figs. 15A and 15B are flow diagrams indicating the steps in the process of 
calculating the correlation value. 

Fig. 16 is a flow diagram indicating the steps in the process of calculating the bound 
15 correlation value. 

Fig. 17 is a flow diagram indicating the steps in the process of reconstructing the 
fullpass image. 

DETAILED DESCRIPTION 

20 Fig. 1 is a schematic illustration of a representative system 10 for both capturing a 

visible image and an infrared image from the frame 12 of film 14 and storing the data in a 
suitable computer 16 for additional processing as further described hereinafter. An 
appropriate sensor, such as a color camera 18, receives image information from the frame 
12 when it is illuminated by a light source, such as visible light source 20 or infrared light 

25 source 22. The process of generating sequential images may of course be automated as 
shown generally by the control line 24 from the computer 16 to the image generating 
components of the system 10 as is well known in the art. A keyboard 26 provides a 
conventional user interface to the computer 16, and a color monitor 28 is further provided 
for viewing the various images as desired. 

30 When implemented, the present invention, unlike the prior art, nulls defects even in 

the case where there are discrepancies between the records. The present invention nulls, 
that is to say eliminates by subtraction, the effects of surface and near surface defects on the 
image media. The processing computer 16 is, in effect, given not only a visible record 
containing both image and defect information but also an infrared record containing only the 

35 defect information. In a preferred embodiment, the infrared record further consists of a 
gain control to adjust the intensity of the information in the infrared record. In practice, the 
gain control needs to be adjusted over a somewhat limited range. This range may only need 
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to be a few percent, such as between zero and five percent. The computer also has the 
ability to adjust the gain conirol to achieve the best nulling effect after the division process, 
such that none of the defect appears in the image, either as a positive ghost or as a negative 
ghost. 

5 An overview of how the gain is used to control the intensity of the information 

contained in the infrared record will be described next. This description will cover use of 
the gain in two different ways. It may be either multiplied by the image information, or the 
image information may be raised to the power of the gain. 

In the first type of use, assume that the gain is a constant, which when multiplied by 
10 the information in the infrared record, gives an infrared product. The information in the 
visible record is then divided by the infrared product. The resulting corrected image will be 
made brighter or dimmer by the gain constant depending on whether the constant is less 
than or greater than one respectively. However, using the gain constant as a simple 
multiplier will produce no resulting effect on the degree of defect nulling since the entire 
1 5 corrected image will have been subjected to the same gain constant. 

In the second type of use, the infrared record data is raised to the power of the gain, 
which produces an indication of "intensity." In the context of nulling, a definition of 
intensity is the contrast measured (that is to say, the ratio of light to dark or gamma). The 
infrared record is raised to the power of the gain. The image data is then divided by the 
20 infrared data after it has been raised to the power of the gain. For example, supposing a 
particular defect absorbs 20% of the light. Without the defect, the visible signal might 
measure 50% and the infrared signal 100%. With the defect the visible signal would 
measure 50%* (100% - 20%) = 40% and the infrared signal would measure 100%* (100% 
- 20%) = 80%. To further facilitate the manipulation of the image data, the visible and 
25 infrared records are converted to natural logarithms. Use of logarithms is expedient 
because it not only changes the division of linear visible data by the linear infrared data to a 
subtraction of the log infrared data from the log visible data, but raising the linear infrared 
data to the power of the gain, or gamma, becomes a multiplication of the gain by the log 
infrared data. 

30 As a further refinement, the gain may be adjusted separately for different spatial 

frequency bands. This overcomes- some of the limitations in real lenses and CCD sensors 
that receive an infrared image with less focus and sharpness than a visible image. For more 
effective nulling of defects, the computer is programmed to find a higher gain matching the 
infrared record at higher frequencies to overcome this limit in real lenses and sensors. 

35 Without this different correction for spatial frequency, an effect is seen similar to that 
produced by an unsharp mask in which the details of a defect are attenuated less than the 
overall defect. Further, it has been found that the degree of correction required tends to 
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vary between scans as the focus errors in the lens interact with film warpage. Hence, 
nulling must be performed on each specific image, and even on specific areas in an image, 
rather than as a rote sharpening of the entire infrared image. 

A cross correlation or autocorrelation effectively averages the power of all 
5 frequencies; however, if the mean has a non-zero magnitude, this zero frequency mean, or 
"DC" term, usually predominates. Because the mean is not part of the signal of interest in 
this process, the cross correlations and autocorrelations would give false indications. The 
division of the visible and infrared images into frequency bands has an advantage in that 
each band has zero mean, solving the potential problem of the zero frequency mean with the 
10 autocorrelation and cross correlation. However, even without a nonzero mean, it is still 
possible for widely different frequencies to beat together in the cross correlation and 
autocorrelation terms similar to heterodynes in radio theory. Use of discrete frequency 
bands has the additional advantage of preventing widely different frequencies from 
combining in this way. Frequency bands also allow the use of the logarithmic domain to 
15 simplify the application of a gain and the division of visible pixel data by infrared pixel data 
as a subtraction while still allowing the cross correlations and autocorrelations to be 
performed with zero mean. 

In another embodiment, gain is adjusted separately for different spatial regions of 
the image. This has been found to overcome a number of limitations in real scanners that 
20 illuminate and focus the edges of a film differently than the center. Further, use of spatial 
regions overcomes problems with different kinds of defects exhibiting slightly different 
attenuation characteristics between visible and infrared light, often because of a 
misalignment between the scanner's visible and infrared illumination systems. The region 
over which the gain is allowed to vary is called the "smudging" region, 
25 The computer is programmed to calculate nulling gain in one of several ways. In 

one embodiment, the computer generates a measure of "likeness" over a smudging region 
by measuring the cross correlation between a candidate restored image and the defect 
image. The cross correlation is the product of the image and defect, pixel by pixel, called 
the point cross correlation, averaged over the smudging region. If the cross correlation is 
30 zero, then that defect has been nulled. If the smudging region is too small, the process 
might null out a portion of the image that looks like the defect along with the defect; 
however, if the smudging region is large enough, the chance of image detail looking 
holographically like the defect is remote, and the image is preserved as the defect is nulled. 
In another embodiment, the nulling gain is calculated directly using the cross 
35 correlation between the uncorrected visible image divided by the autocorrelation of the 
infrared image. First, the cross correlation is averaged over the smudging region. Then the 
autocorrelation, which is the square of the value of each pixel (called the point 
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autocorrelation), is averaged over the smudging region. The average cross correlation is 
divided by the average autocorrelation to produce an autocorrelation gain for the region. 
When this autocorrelation gain is applied to the infrared image and canceled from the 
visible, a corrected image is formed that gives a zero cross correlation with the infrared 

5 image. , 

The process of the preferred embodiment will now be summarized. First, both 
linear visible and infrared scans of the image are taken in a conventional manner as 
described in U.S. Patent No. 5,266,805. Pixel data from both scans are converted into 
logarithmic form, and then each of these logarithmic images are split into frequency bands. 

10 For every corresponding frequency band, the pixel by pixel cross correlation between the 
visible and infrared records is calculated as is the autocorrelation of the infrared pixel data. 
Next, the pixel by pixel cross correlations and autocorrelations are averaged, or 
"smudged," over small regions. For each region, the smudged cross correlation is divided 
by the smudged autocorrelation to find the nulling gain. The log of the infrared bandpass 

15 data is then multiplied by the nulling gain, and the product is subtracted from the 
logarithmic visible image data to obtain the corrected image data for that bandpass. After 
this operation has been performed on all regions of the image and all bandpasses, the 
corrected bandpasses are reassembled to obtain the logarithmic data for the corrected image. 
The corrected image may then be displayed, stored, or transmitted after the antilog of the 

20 data is calculated. 

Referring now to Fig. 2, a generalized description of the overall process of defect 
channel nulling is provided to facilitate its understanding. At the beginning of the process, 
the dimensions of the input image and defect data are ascertained at step 40, and program 
variables are initialized at step 42. Both the image record 44 and defect record 46 are 

25 converted into a logarithmic form by steps 48 and 50, respectively. From the log image 
record 52 and log defect record 54, the log image bandpass pyramid 56 and log defect 
bandpass pyramid 58 are created by steps 60 and 62, respectively. Layers of the image and 
defect bandpass pyramids 56 and 58 are then correlated by step 64 to create a corrected log 
image pyramid 66. Based upon the correlation, the image pyramid 66 is corrected and 

30 reconstructed at step 68 to a fullpass image 70. Finally, the corrected log image 70 is 
converted by step 72 from the logarithmic format 70 to a corrected linear format 74. Each 
of these high level steps will now be described in detail below. 

Fig. 3 is the first of the diagrams which detail the processes more generally outlined 
in Fig. 2. Fig. 3 provides a detailed description of the variable initialization of the pyramids 

35 mentioned in Fig.2 to begin the overall defect channel nulling process of the present 
invention. First, the sequence of the array is set at 0 by step 80 to begin the counting 
process. Next, the image width and height are defined by step 82. As the levels of the 
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array are increased, the width and height of the image are calculated as a function of the 
levels in the array at step 84. This process is repeated as long as the width and height of the 
current level are not both equal to one as checked by condition 86. When the width and 
height of the image are both equivalent to one at decision step .88 thereby indicating a 

5 pyramid level consisting of a single pixel, the counting of the number of levels is completed 
and is set by step 90. With this step, variable initialization is complete. 

In one embodiment, the next step is to convert the image and defect records into 
logarithmic format. Converting to a log format is desirable because of the pyramid building 
process which facilitates splitting the images into frequency bands. Each layer of the 

10 pyramid represents a frequency band. Using the pyramid process allows a zero expected 
value, or mean, at each level and therefore eliminates the bias or "DC" term described 
earlier that differentiates a covariance and correlation. Because a zero mean can be reached, 
a linear format which could possibly result in a division of zero by zero (if the visible 
record is differentiated from the defect record by a division) cannot be used. Converting to 

15 a logarithmic scale solves this problem because division of the image record by the defect 
record is effected by a simple subtraction between equivalent levels of the image and defect 
pyramids. Thus, use of the log domain combined with the pyramid process solves 
problems both with DC bias in the correlation calculations and with effective division on 
separate frequency bands. 

20 When using the disclosed method to null out defects, it is helpful to be able to 

measure how the spatial frequency response of the defects in one record compares to the 
response in another record. One way to do this is to perform a fast Fourier transform on 
each record and compare them. Based on the comparison, the attenuated frequencies of one 
record can be boosted to match the other. Another expedient is to perform two-dimensional 

25 cross correlation and autocorrelation calculations. By spatially dividing one by the other, 
relative blurring in the image may be found and corrected. In a preferred embodiment, each 
image is divided into frequency bands. By performing the nulling process on each 
frequency band separately, any variations in frequency response between the images will be 
automatically compensated thereby boosting any attenuated frequencies. 

30 One way of dividing the frequency into frequency bands that allows the original 

image to be reassembled is to create a band pass pyramid. In a bandpass pyramid, the 
image is first low pass filtered and then subsampled at a lower resolution. The lower 
resolution is then upsampled to the original resolution in order to produce an image devoid 
of details that could not be represented by the lower sampling frequencies but which 

35 contains all the details that could be represented by the lower frequencies. When this 
upsampled image is subtracted from the original image, only the high frequency details 
remain. The entire process is repeated on the subsampled image just discussed and the 
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process is continued recursively to generate successively lower and lower frequency bands. 
Each band represents only those details that could not be represented by the lower 
frequency band. Because of the way that each of these bands is generated, if the lowest 
frequency band is upsized and added to the band above and the procedure is repeated 
5 recursively until the highest band is reached then the original image will be precisely 
reconstructed. Possible variations in this process include the way in which an image is 
downsized or upsized. If a poor algorithm is used, there will be aliases in the bandpasses. 
These aliases will cancel each other out as long as nothing is done to each band; however, if 
the different bands are processed with different gains, the aliases will affect the final 
10 reconstructed image. Thus, the upsize and downsize chosen in a preferred embodiment are 
a compromise between minimizing aliases and conducting fast computations. Other 
downsizes and upsizes which are possible include the bicubic process used in Photoshop 
manufactured by Adobe Corp. 

Fig. 4 describes the pyramid building process in greater detail. Typically, a 2,000 ( 
15 3,000 original pixel image is used and is downsized by one-half at each level. For 
illustrative purposes, an original 8(8 image 100 is used. In the first step, the image 100 is 
downsized to a 4(4 image 102. Figs. 5A through 5D provide a detailed description of the 
process of downsizing. They show how a 3(3 region 104 at one level of the original 8(8 
image 100 is downsized to a single pixel 51 at the next level. Each pixel in the 3(3 region 
20 104 is multiplied by a factor which is based upon the pixel's position within the 3(3 region 
104. In Fig. 5A, the pixel with designation LI is given a factor of one, pixel L2 is given a 
factor of two, L3 has a factor of one, L6 has a factor of two, L7 has a factor of four, L8 
has a factor of two, LI 1 has a factor of one, L12 has a factor of two and L13 has a factor of 
one. Each pixel and its corresponding factor are shown in diagram 106. These factors are 
25 chosen so that the resulting single pixel is a weighted average of the pixels in the other level 
based on the proximity of the source level pixels closest to the resulting pixel. Next, the 
sum of the product of pixel values and factors is taken and divided by sixteen, which is the 
sum of the factors so that the value for the resulting pixel is anti-aliased. The value after the 
division is the value of pixel S 1 which is the downsized representation of the original 3(3 
30 region 104. Use of the factors followed by the subsequent division by the sum of the 
factors permits each pixel in the upsized level to contribute a weighted portion to the value 
of the pixel being calculated in the downsized level based on the upsized pixel's proximity 

to the downsized pixel. 

The remaining diagrams in Figs. 5B through 5D are examples of the above 
35 described process over different areas of the original 8(8 image 100. Fig. 5B illustrates the 
same 8(8 area 100 with a different 3(3 region 108 within the 8(8 area 100 being sampled. 
As in Fig. 5A, the pixels of Fig. 5B are multiplied by a factor based on each pixel's 
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respective position in the 3(3 region 108. The pixel with designation L3 has a factor of 
one, L4 has a factor of two, L5 has a factor of one, L8 has a factor of two, L9 has a factor 
of four, L10 has a factor of two, L13 has a factor of one, L14 has a factor of two, and L15 
has a factor of one. Again, the sum of the pixel-factor products is taken and divided by 

5 sixteen, the sum of the factors. The value after the division is the value of pixel S2 which 
is the downsized representation of the 3(3 area 108. Fig. 5C illustrates the same 8(8 area 
100 of Figs. 5A and 5B with a different 3(3 region 1 10 being sampled. As in Figs. 5A and 
5B, the pixels of Fig. 5C are multiplied by a factor based on each pixel's position in the 3(3 
region 110. The pixel with designation Lll has a factor of one, L12 has a factor of two, 

10 L13 has a factor of one, L16 has a factor of two, L17 has a factor of four, L18 has a factor 
of two, L21 has a factor of one, L22 has a factor of two, and L23 has a factor of one. The 
sum of the pixel values is taken and divided by sixteen which is the sum of the factors. The 
value after the division is the value of pixel S4 which is the downsized representation of the 
3(3 area 1 10. Fig. 5D illustrates the same process with a different 3(3 region 1 12 of the 

15 8(8 area 100 being sampled. The value after the factor multiplication, summation and 
division calculations is the value of pixel S5 which is the downsized representation of the 
3(3 area. 

Returning to Fig. 4, the new downsized image 102 is then upsized to create an 
upsized image. Fig. 6 illustrates a region of a 4(4 area being upsized to an 8(8 area 1 18. 

20 There are several ways to upsize including using a sine function (sin x divided by x). 
However, the algorithm chosen for this invention addresses the calculation of a higher 
resolution pixel by looking at the pixels in the downsized image lying in four directions 
from the higher resolution pixel. This algorithm also has the advantage of being more 
expedient in that it is much faster to compute than a sine function. 

25 The pyramid levels may be conceptually considered to lay over each other. Each 

covers the same image, but a lower level covers it with pixels more sparsely. Pixels in a 
specific level have pixels both over and under them in higher and lower levels. In addition, 
as seen in Fig. 7A through 7D, pixels also have pixels to the left, to the right, to the top, 
and to the bottom in the same and in different levels. As used herein, "over and under" will 

30 refer to the relation between different levels while "to the top" and "to the bottom" will refer 
to a spatial relationship within the same image level. In the case where the value for the 
pixel to be calculated at a higher resolution has a pixel directly under it at a lower resolution, 
the lower resolution pixel is simply copied to give the higher resolution pixel. A different 
situation occurs when the higher resolution pixel being calculated has lower resolution 

35 pixels under it located to the right and to the left in the lower resolution level but not directly 
under the higher resolution. In this latter case, the algorithm averages the two pixels of the 
lower level which lie to the right and left under the high resolution pixel. A third situation 
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arises when the higher resolution pixel to be calculate* has lower resolution pixels £*r * 
, , t0 the top and to the bottom in the lower resolution level, but no, dtrectlv under «. In 

case, Z aigorithm avera g es the two .ower resolution pixels to £ve 
resolution pixel. A fourth scenario is where the higher resolution ptxel bemg calculated 
; L not have any pixe!s under it, whether to the right, to the left, to the top or to ft. 
bTrL, or direcuy under it. Here, the higher resolution pixel is directly .» between the 
lower resolution pixels and the algorithm takes an average of four pixels. 

Hgs. 7A ftrough 7D describe upsizing in general and il.us.ate the four cases 
described Fig. 7A focuses on Pixel SI. Pixel SI is directly upsized to post-on U7 a. 120 
0 "gthedatainS.toU7. Pixels SI and S2 are summed and divided by two, ogtve 
Le U8 equal to me average of SI and S2 at 122. Snarly, p.xe U SI and S4 ^ 
led and divided by two to give value U12 a, .24. 

summed and divided by four to give the resulting average to ptxel va ue U13 at 126^ to 
; g. 7B, the upsizing is focused on pixel S2. Pixel S2 is upsized direct,, to p«- W« 
5 ,28. Pixels S2 and S3 are summed and divided by two to g,e value U10 a, 130. ftxe 
S2 and S5 are summed and divided b, two to give value U14 a, ^132 Ptx Is S2 S3. S5 
m d S6 are summed and divide, by four ,0 give pixel value U.5 a. ,34 In F* TC. «. 
upsizing is focused on pixel S4. Pixel S4 is upsized directly to posmonU! at IX. 

are summed and divided by four to g.ve pixel value U23 at 142 InMg , " 
focused on pixel S5 and upsizes to a different region. ValuesforU19,U20,U24andU25 
arecalculatedexactiyasdescribedarx W einconnectionw.d.F 1 g.7C. 

Referring back to Fig. 4, the new upsized image (not shown), wh.ch came from 

bandpass image .44. F.g. 8 Ulustrates in detail the different process. F*ch e» 
fte upsized smaller image .48 is subbed from the coupon drf o f « 
image 150 to produce an element in bandpass image 146. For example, BP1 equals U 
mTs U7, J BP2 e,ua,s L2 minus US. As shown in Fig. 4 this process » , rep«£ 
,0 recursively such mat each btmdpass image represents a specif, frequency ta d of ft. 
origina. image 100. The lowest level of the pyramid » a stngle p.xel 152 whtch ts 
bandpassed because there is nothing from which to difference,,. . ^ Dlssed 

Ate the original image has been bandpassed, the defect ~ 
van. the same pyramid building process as described above. As wtth tire image, each level 
35 Ihltf^ZMrepresetttsafreauencybandofuteoriginaldefect. ^etmage 

associated image level to recover an image without defects at that level. Pnor to subtraction 
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from the corresponding image level, the defect level is multiplied by a gain chosen to match 
the defect part of the image to provide a better null after subtraction. Preferably, this gain is 
made a function of position within the defect record so that each element of the defect 
record can be adjusted with a different gain. 
5 There are many ways of finding this gain. One way is to iteratively try different 

gains until detail in the subtracted image record does not match detail in the defect record. 
In other words, the gain is chosen based on which one makes the image record and defect 
record as uncorrected as possible. The defect channel nulling of the present invention uses 
a non-iterative method involving direct calculation to reach an adequate degree of detail 
10 uncorrelation. One measurement of the detail matching is the mathematical crosscorrelation 
in which the two records are multiplied element by element. To the degree both records 
match, their detail will track each other closely, in either case producing a positive result 
after the element multiplication. To the extent that they are uncorrected, the 
crosscorrelation multiplication will sometimes produce a positive product, but will just as 
15 likely produce a negative product. Thus, when the crosscorrelation is averaged over a 
region, any nonzero value indicates a match or correlation. 

Other refinements in gain are based on prior knowledge of record characteristics. 
- For example, in one embodiment, the image record is set to zero whenever the defect rises 
above a predetermined intensity. Setting the image record to zero when the defect intensity 
20 is high is an acknowledgement that beyond a certain point, an accurate prediction based on 
the image is impossible, and the best estimate is just zero. Although it is possible to set the 
image record pyramid level to zero whenever the magnitude of the defect record exceeds a 
threshold, the present embodiment uses instead a refinement that allows a gradual zeroing 
or blending so that no abrupt drops to zero are created in the image record. Another 
25 possible refinement approach includes varying the size of the averaging area based on 
image and defect characteristics. Thus, the gain refinement used by the present 
embodiment is given as an example, and is not intended to limit the process to a specific 
mode of calculating the gain. 

The next step in the defect nulling process is to correlate the bandpass frequency 
30 band for both the image and the defect pyramids. An overall description of the correlation 
process is illustrated in Figs. 9A and 9B. The correlation process begins with the 
generation of a correlation histogram at step 156 for each corresponding level in the image 
and defect pyramids. Figs. 9A and 9B give a detailed description of how to generate a 
correlation histogram for a single level of the image and defect pyramids. Processing starts 
35 with each image pixel and the corresponding defect pixel at a specified level. First, 
variables used in the process are initialized at step 158. Specifically, the inner limit for the 
histogram is set to zero and the outer limit is set to two. These limits are chosen to center 
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around the ideal situation when the correlation equals one. Because the correlation is a 
squared number, its lowest value is zero, thereby setting the inner limit for the histogram. 
In order for the approach to be balanced around one, the upper limit is set at two. 
However, it is understood that other limits could be used. 

5 After the histogram is generated, the next step 160 is to calculate the correlation 

cutoffs which are the new upper and lower cutoff limits for the histogram. 

A graph of the histogram with the correlation cutoffs is illustrated in Fig. 1 1. This 
figure shows how well the defect record tracks the visible component of a defect. The 
x-axis 162 is the ratio of the defect record to the visible record when the visible image has a 

10 defect. The y-axis 164 is the number of pixels representing that ratio weighted by the 
square of the pixels. If, for example, the visible defect exactly tracked the defect record, all 
pixels would show a one to one correspondence. In practice, however, the infrared and 
visible records do not track exactly, and the graph of Fig. 11 shows the degree to which 
they diverge. The ten percent points 166 and 168 as illustrated attempt to show the range 

15 for which articulation or fine adjustment of defect gain relative to visible is expected to be 
needed. The process of calculating the correlation cutoffs is described in Figs. 12A and 
12B. These new limits 166 and 168 define the range within which the values will be 
considered. 

After the cutoff values are defined, a crosscorrelation and autocorrelation are 

20 calculated as illustrated in Fig. 13. Applying the cutoff values from Fig. 12B to the image 
(also referred to as preclamping) produces a clamped image in which each individual image 
element is limited to a present value. Fig. 14 illustrates the process of preclamping which is 
a way of expressing the importance of a particular pixel. Normally, it is expected that a 
defect in the visible record is about 1.0 times the defect record. The purpose of nulling is to 

25 accommodate those cases that deviate from 1.0. The present embodiment preclamps the 
visible record to between 0.0 and 2.0 times the defect record prior to taking the 
crosscorrelation and autocorrelation. In this way, the effect of large variations in the visible 
record are taken into account but are limited in value before the multiplication of the 
crosscorrelation and autocorrelation are performed. In Fig. 14, this is done by setting the 

30 variables Low Clamp equal to 0.0 and Upper Clamp equal to 2.0. If the visible record 
were, for example, 4.0 times the defect record, the preclamping would treat the situation as 
though it were only 2.0 times the defect record, but would maintain full weighting, which 
is mathematically equivalent to letting it be 4.0 while only counting it with half the weight. 
Another weight which can be used to express confidence in a pixel is based on how far the 

35 calculated correction coefficient (defined below) varies from the expected value of unity for 
each individual pixel. If the variation from unity is far, that pixel is given little weight and 
vice versa. In a specific example, any pixel that gives a single pixel correction coefficient 
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between 0.0 and 2.0 is given full weight. For a correction coefficient beyond these values, 
the preclamping process effectively reduces weighting of the particular pixel proportional to 
just how far the correction coefficient is beyond these limits. Upon completion of the 
process described in Fig. 14, the clamped element rather than the element itself is then used 
5 to calculate the crosscorrelation and autocorrelation. 

As shown in Fig. 13, the crosscorrelation between the image and the defect record 
is, for a single pixel, the product of two values, the clamped element multiplied by the 
defect element. The defect autocorrelation for the same pixel is the defect element value 
squared. If the crosscorrelation is then divided by the autocorrelation for that single pixel, 
10 the result, in theory, is the correction coefficient which is the amount of the defect record 
present in the image record. However, if this approach is strictly applied, the entire image 
will be nulled to zero for each pixel, which is obviously undesirable. Instead, only that 
component of the image that matches the defect should be nulled. This is done by 
averaging these values over a region. This averaging process or "smudging" is described 
15 in Fig. 9B. In one embodiment of the invention, a 3(3 element area is used for the 
averaging region because it encompasses an average full wave in the pyramid frequency. A 
larger area would give a more accurate estimate of the correction coefficient; however, 
. smaller areas allow better articulation around the edges of defects which, because of 
imperfections in the angle of lighting between different colors, may require the correction 
20 coefficient to vary pixel by pixel. 

After the single element cross and autocorrelations for each element of the 3(3 
element region are calculated beginning at step 190, the nine cross correlations are averaged 
and the nine autocorrelations are averaged at step 192. As previously mentioned, this 
averaging process across a region is called smudging. The average of the cross correlations 
25 is termed the regional crosscorrelation value. The average of the autocorrelations is the 
regional autocorrelation value. Finally, the regional crosscorrelation divided by the regional 
autocorrelation to give a gain measuring correlation. For an individual pixel approaching 
zero defect where both the autocorrelation and crosscorrelation are zero, little is contributed 
by the pixel to the overall average for its region. On the other hand, at a pixel with a strong 
30 defect record, both the crosscorrelation and autocorrelation are large, and the corresponding 
pixel will contribute significantly to the average for its region. 

From the regional cross correlation and regional autocorrelation values, the 
correlation value is determined at step 194 in Fig. 9B. This process is further detailed in 
Figs. 15A and 15B. First, an articulation gain value which is the exact gain at a pixel 
35 required to null the visible record is calculated at step 196 by dividing the image element by 
the defect element. This is equivalent to dividing the single element crosscorrelation, or 
image element multiplied by the defect element, by the single element autocorrelation, or 
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defect elemen. squared. As previously described, the correlation value is obained by 
dividing the regional crosscorrelation by .he regional autocorrelation a. step 198. The 
threshold mask value is obtained by dividing .he regional au,o«>rrelauon by a 
predetermined threshold constant a. step 200. In this embodiment, the threshold constant ,s 

5 015 U^m B om^^«^^«^'>^^^^ ta ^ P °r 
then that pixel is allowed to "articulate" the required correction coefficient for complete 
auction the visible record a. that particular pixel, and change the visible tmage pixel 
by pixel, wi thout regional averaging. Note that if the threshold mask is zero a, step 202, 
tine new upper cutoff will be equal to the new lower cutoff which is simply equal to the 

10 correlation as shown a. step 204. If the threshold mask is one a. step 206, the new upper 
cutoff equals the old upper cutoff, originally assumed to be 2.0. and the new lower cutoff 
equals .he old lower cutoff, originally assumed to be 0.0 as seen a. step 204. Correlation .s 
set to articulation a. step 208 if articulation is between new upper cutoff and new lower 
cutoff, but is limited to the range between new upper cutoff and new 'o^utoff In 

15 Act. if the defect intensity is low, correlation will be left unchanged, but tf the defec 
intensity is high, then chelation will be equal to articuiation. In other words, ,f*e defec 
intensity is high, ten the correlation will track articulation for complete nulhng of the unage 

in the visible record. _ _ 

The range of correction is still limited, however, to be close to the average to 
20 preventexrremecancellationsoutsidetherangeofprobability. Using me threshold mask 
value, new upper and lower cutoff values are generated which serve to bound the 
correlation value generated above, mis process is referred to in Fig. 9B . ;ca— £ 
the bound correlation value a, step 210 and is described in deuttl » Ftg. 16. The upper 
cutoff arai lower cutoff values may be made more restrictive than the full range 0.0 to 2.0 
25 in particular, it has been found to improve reconstruction quality when only the largest 
pyramid levels have the full range, and the smaller levels are constrained » a narrow range 
based on values of correlation in corresponding areas of higher levels. F.g . 1 i iU*~ s 
an algorithm to set the upper cutoff and lower cutoff to be plus and minus 20% from the 
actual correlation found for the previous larger pyramid level. 
30 As shown in Fig. 9B, this new bound correlation value is ten used to correct the 

image element at step 2!2. The new correlation value is multiplied by the defec, element 
and mis product is subtracted from the image element Tne image element ts now fully 
corrected. After each image elemen. is corrected in the manner described, re— ..on of 
the fullpass image begins as illustrated in Fig. .7. Starting from the bonom of the ur^ 
35 pyramid, each level is upsixed by adding each elemen. of the upazed level to each element 
eT fa level below. Ate all .he levels have been added, me reconsttuction of ft. log unage 
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is complete. Upon completion, the antilog is taken of the log image and the resulting image 
can now be displayed, stored, or transmitted. 

While this invention has been shown and described with reference to particular 
embodiments thereof, it will be understood by those skilled in the art that the foregoing and 
5 other changes in form and detail may be made therein without departing from the spirit and 
scope of the invention. 
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9. The method of claim 1 wherein the defect signal includes a plurality of defect signal 
frequency bands, and wherein the gain is adjusted differently between at least two of the 
bands. 

5 10. The method of claim 9 wherein one of the defect signal frequency bands comprises 
the D = Factor of the image signal. 

1 1 . The method of claim 9 wherein the defect signal frequency bands are logarithmically 
spaced by a frequency incrementing factor. 

10 

12. The method of claim 1 1 wherein the frequency incrementing factor is two. 

13. The method of claim 9 wherein the defect signal frequency bands are linearly 
spaced. 

15 

1 4. The method of claim 9 wherein the step of adjusting the gain to minimize the defects 
in the recovered image signal comprises the step of calculating the correlation between the 
image signal and at least one of the defect signal frequency bands. 

20 15. The method of claim 14 wherein the step of adjusting the gain further comprises the 
step of calculating the gain as a function of the correlation between the image signal and at 
least one of the defect signal frequency bands. 

1 6. The method of claim 9 wherein the defect signal received is further characterized to 
25 include a plurality of defect signal regions, and wherein the gain is adjusted differently 

between at least two of the defect signal regions. 

17. The method of claim 16 wherein the step of adjusting the gain to minimize the 
defects in the recovered image signal comprises the step of calculating the correlation 

30 between the image signal and at least two of the defect signal regions. 

18. The method of claim 17 wherein the step of adjusting the gain comprises the step of 
calculating the gain as a function of the correlation between the image signal and at least two 
of the defect signal regions. 

35 

1 9 . The method of claim 1 8 wherein the step of calculating the gain for a selected defect 
signal region within a selected defect signal frequency band comprises the steps of 
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means for receiving from the image storage medium a defect signal, in register with 
the image signal, indicating the defects in the image; 

means for applying a gain to the defect signal to generate a normalized defect signal; 
means for differencing the image signal and normalized defect signal to generate a 
5 recovered image signal; and 

means for adjusting the gain to minimize the defects in the recovered image signal. 

26. The apparatus of claim 25 wherein the recording medium is a substrate bearing a 
two dimensional image. 

10 

27. The apparatus of claim 26 wherein the substrate contains a visible image. 

28. The apparatus of claim 26 wherein the defect signal is received from the substrate in 
infrared. 

15 

29. The apparatus of claim 25 wherein the gain adjusting means includes means for 
measuring the defects in the recovered image signal. 

30. The apparatus of claim 29 wherein the measuring means comprises means for 
20 calculating the correlation between the recovered image signal and the defect signal. 

31. The apparatus of claim 25 wherein the gain adjusting means comprises means for 
altering the gain as a function of the correlation between the recovered image signal and the 
defect signal. 

25 

32. The apparatus of claim 25 wherein the gain adjusting means comprises means for 
measuring the defects in the image signal. 

33. The apparatus of claim 32 wherein the measuring means comprises means for 
30 calculating the correlation between the image signal and the defect signal. 

34. The apparatus of claim 33 wherein the gain adjusting means comprises means for 
calculating the gain as a function of the correlation between the image signal and the defect 
signal. 

35 
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f i„, m 35 wherein one of the defect signal frequency bands is a 
5 36. The apparatus of claim 35 wnerein u« 

mean. 

« ■ m « wherein the defect signal frequency bands are 
37 The apparatus of claim 35 wherein in 
togarithmically spaced by a frequency incrementing factor. 

spaced. 

^eans for calculating the correlation between the image sign 
signal frequency bands. 

. at least one of the defect signal frequency bands. 

* u •« the defect signal is further characterized by 
42. The apparatus of claim 35 wherein the drf*. g 
25 comprisingapluralityofdefectsignalre^andwheremth 

between at least two of the regions. 

30 registered with the visual image data, compns^ ^ ^ ^ 
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respective spatial frequency bands; non -image defect data sets for 
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35 the visual imaging data; 
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cross correlating the visual image data set and the non-image defect data sets for 
respective spatial frequency bands to obtain crosscorrelation variables for each of the spatial 
frequency bands; 

autocorrelating the non-image defect data set for each of the spatial frequency bands 
5 to obtain autocorrelation variables for each of the spatial frequency bands; 

. processing the crosscorrelation variables over a plurality of spatial regions to obtain 
spatially defined crosscorrelation variables for each of the spatial frequency bands; 

processing the autocorrelation variables over a plurality of spatial regions 
substantially similar to the spatial regions of the crosscorrelation variables to obtain spatially 
1 0 defined autocorrelation variables for each of the spatial frequency bands; 

processing the spatially defined crosscorrelation variables with the spatially defined 
autocorrelation variables to obtain spatially defined nulling gains for each of the spatial 
frequency bands; 

applying the spatially defined nulling gains to the respective non-image defect data 
1 5 sets to obtain nulling data sets for the respective spatial frequency bands; 

applying the nulling data sets to the visual image data sets for the respective spatial 
frequency bands to obtain corrected visual image data sets for the respective spatial 
frequency bands; and 

merging the corrected visual image data sets for the respective spatial frequency 
20 bands to obtain corrected visual image data. 

44. A method of processing visual image data with non-image defect data to remove 
non-image defect information from the visual image data, the non-image defect data being 
registered with the visual image data, comprising: 

25 comparing the visual image data with the non-image defect data to determine a 

nulling gain; 

applying the nulling gain to the non-image defect data to create nulling data; and 
applying the nulling data to the visual image data to remove non-image defect 
information therefrom. 

30 

45. The method of claim 43 further comprising identifying a plurality of spatial 
frequency bands to overcome limitations in real optics and image sensors, the visual image 
data further comprising a plurality of respective visual image data sets for the spatial 
frequency bands, the non-image defect data further comprising a plurality of respective 

35 non-image defect data sets for the spatial frequency bands, the nulling gain further 
comprising a plurality of respective nulling gain sets for the spatial frequency bands, and 
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<LowerCutoff index - 
LawSuift = Yhist 



= YHISTSIZE - 
[LowState] 



> 




LmvState = LowStote - 1 
LowSum = LowSum + Yhist [LowState] 



Lowe rCu toft index = LowState 
LowerCutoff = Range * (LowState / YHISTSIZE) 



FIG. 1ZB 
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Calculate CrossCor relation and Autocorrelation 



(Calculations A 
begin hereby 



Clamped Element = Calculate clamped element 



CrossCorrelctior. » C taliped Element * defect 
element 

Autocorrelation = defect element * defect element 



FIG. 13 



WO 98/31142 



Calculate Clamped Element 
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Receive: 
image element 
defect element 
low c lamp 
upper clamp 



( Clamping ^ 
V^begins hare./ 



Clamped element = image 
element 



No/ defect ^^Yes 
element > 0 




ClonpedEleneiil = | ( 

LawClai^deffict olooent 



Ciaipediraage Ss, ^Yj s 
AJpperC lamp 



ClifflpedElanient" 
LowC lamp *de feet element 



CI(npedElBflBnl = 
UpperCluip*defect eleeeni 




ClaipedElenent = 
UpperClanp^defect ekient 



FIG. 14 



Done. 
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Calculate Correlation value 



196-x 



Correlation calculation 
3tarts here. 



Articulation = image element / defect element 



Correlation = Reg ionalCrossCor relation/ 
Reg 10 na lAutoCo r re la t ion 



200 



ThresholdMask = ReglonalAutoCorrelation/ 
THRESHOLD 



204 



1 




202^ 




Thresha IdWask = 0 



ItevlTppafCutaff - (lipparCutdfrthrBstoldtosk) + (Correlation t (1 - ThreshofdMosk)) 
NewLowerCutoff - (LowerCu lof f*Th r«sho ItUnsk) f (Correlation ♦ (1 -Thresho IdMosk)) 



FIG. 15A 



'3/25 




FIG. 15B 
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Bound Correlation Value 
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Bounding of 
Correlation 
^begins hereb y 



AvgPrevCorr = Average of previous level Correlation values 



UpperAvgPrevCorr = AvgPrevCorr + (20% of AvgPrevCorr) 
Lower AvgPrevCorr = AvgPrevCorr + (20% of AvgPrevCorr) 



Yes 




UpperAvgPrevCorr = UpperCutoff 



Yes 



LowerAvgPrevCorr < LoiwrCutoff 



LowerAvgPrevCorr = LowerCutoff 
I 



Yes 




Correlation - UpperAvgPrevCorr 



Yes 



Correlation = LowerAvgPrevCorr 




Done. 



FIG. 16 



WO 98/31142 



fCWlTS97/24i36 



25/25 



Reconstruct Fiillpasa Imcqe 



I Reconstruction j 
^ begins here, J 



■ ^i = NiimOfLevels-l^ 



For each level of the image pyramid starting from 
the 'bottom' (i.e. the level with dimensions 1x1). 



I 



Upsize level[i] 



Add each element of the up si zed 
levelfi] to each element of level[i-1] 
end save the result in level[i-1] 



i - i-1 




Yes 



Done, The final 
reconstructed image 
is level[D] 
image pyramid 



FIG. 17 
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